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ABSTRACT 



"gj I Two models of the gamma ray burst population, one with a standard candle 

pLn I luminosity and one with a power law luminosity distribution, are x^-fitted to the union 

0> ' of two data sets: the differential number versus peak flux distribution of BATSE's 

long duration bursts, and the time dilation and energy shifting versus peak flux 
information of pulse duration time dilation factors, interpulse duration time dilation 
factors, and peak energy shifting factors. The differential peak flux distribution is 



lO . corrected for threshold effects at low peak fluxes and at short burst durations, and 



the pulse duration time dilation factors are also corrected for energy stretching and 



^^ I similar effects. Within an Einstein-de Sitter cosmology, we place strong bounds on 



the evolution of the bursts, and these bounds are incompatible with a homogeneous 
population, assuming a power law spectrum and no luminosity evolution. Additionally, 



f) i' under the implied conditions of moderate evolution, the 90% width of the observed 

Q . luminosity distribution is shown to be ^ 10^, which is less constrained than others 

j^ i indicate that if the redshifts of BATSE's faintest bursts are to be compatible with that 

which is currently known for galaxies, a standard candle luminosity is unacceptable, 
l^ I and in the case of the power law luminosity distribution, a mean luminosity ^ 10^^ ph 



;h ' s ^ is favored. 



Subject headings: gamma-rays: bursts - cosmology: theory 



1. Introduction 



The angular distribution of the gamma ray burst population has been shown to be highly 
isotropic ( Meegan et al. 1992 ; Briggs et al. 199(^ ). This suggests that the bursts are either located 
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in an extended galactic halo (e.g., Paczyhski 1991) or that they are cosmological in origin (e.g., 



Paczyiiski 1986). Recent measurements of time dilation of burst durations (Norris et al. 1994 



1995; Wijers & Paczyhski 1994; however, see Mitrofanov et al. 1996| ), of pulse durations ( [Norris et 



al. 1996a| ), and of interpulse durations (Davis 1995; Norris et al. 1996b) in the BATSE data, as 



well as measurements of peak energy shifting ( pVIallozzi et al. 1995| ), favor the latter explanation. 

Models, both galactic and cosmological, are typically fitted to the differential peak flux 
distribution of BATSE's long duration (Tgo > 2 s) bursts. Furthermore, this distribution is 
typically truncated at a peak flux of 1 ph cm^^ s^^ to avoid threshold effects. Here, we fit two 
models, one with a standard candle luminosity and one with a power law luminosity distribution, to 
not only BATSE's 3B differential distribution, but also to the pulse duration time dilation factors 
(corrected for energy stretching and similar effects) of Norris et al. (1996a), the interpulse duration 
time dilation factors of Norris et al. (1996b), and the peak energy shifting factors of Mallozzi et 
al. (1995). These three independent sets of measurements are shown to be self-consistent in §4. 
(All three are for long duration bursts only.) Furthermore, via the analysis of Petrosian & Lee 
(1996a), BATSE's differential distribution is extended down to a peak fiux of 0.316 ph cm~^ s~^, 
which corresponds to a trigger efficiency of ~ g o'^ BATSE's 1024 ms timescale. 

Together, the differential distribution and the time dilation and energy shifting factors place 
strong bounds on the evolution of the burst population. These bounds favor moderate evolution 
and are incompatible with homogeneity, assuming only minimal luminosity evolution. This result 
is compatible with the analyses of Fenimore & Bloom (1995), Nemiroff et al. (1996), and Horack, 
Mallozzi, & Koshut (1996). Furthermore, under these conditions of moderate evolution, the 90% 
width of the observed luminosity distribution is shown to be less constrained than others have 
demonstrated it to be assuming no evolution (see §5). Finally, redshift considerations indicate 
that if the redshifts of BATSE's faintest bursts are to be compatible with that which is currently 
associated with the formation of the earliest galaxies, the mean luminosity of the bursts should be 
10^^ ph s~^ or lower. 



r\j 



2. Cosmological Models 

Both the standard candle luminosity model and the power law luminosity distribution model 
assume a power law redshift distribution, given by 

n(z)=no(l + z)^, (1) 

where n{z) is the number density of bursts of redshift z. This distribution is bounded by 

< z < ZM, where zm is the maximum burst redshift. The luminosity distributions of the two 

models are given by 

( (j)o6{L — Lq) (standard candle) 

'^^^ = Uo^-^ (power law) ' ^^' 



The standard candle is of luminosity Lq and the power law luminosity is bounded by minimum 
and maximum luminosities Lm < L < Lm- All luminosities are peak photon number luminosities 
and all fluxes are peak photon number fluxes (measured over BATSE's 50 - 300 keV triggering 
range); however, see recent papers by Bloom, Fenimore, & in 't Zand (1996) and Petrosian & Lee 
(1996b) which introduce the fluence measure. 



2.1. Integral Distribution 

Assuming a power law spectrum and an Einstein-de Sitter cosmology, the bursts' integral 
distribution, i.e. the number of bursts with peak fluxes greater than an arbitrary value F, is given 
for either model by (|Meszaros fc Meszaros 1995| ) 

iV(> F) = ^^ /'" mdL r\l - xf-^-A, (3) 



m 



-^L 



where 



and 



Xo = 'min{xi,X2), (4) 

Xi = r> (5) 

X2 = 1 - ^--T- (6) 

(1 + zm)2 

A photon number spectral index of -1 (or a power-per-decade spectral index of 1) has been 

assumed. This value is typical of burst spectra, especially at those frequencies at which most of 

the photons are received (e.g., Band et al. (1993)). In the case of the standard candle model, eq. 

H becomes 

N{>F)^ / (l-x)'-2Vdx, (7) 

Jo 

where L = Lq in eq. |5[ The factor of proportionality has been dropped because only normalized 

integral distributions (see §3.1) and ratios of integral distributions (see §2.2) are fit to. Eq. has 

the analytic solution 

N{>F)cxf{xo,8-2D), (8) 

where 

^ 2(1 - (1 - x)=^+g) _ 2x(l-x)^+'^ _ XHI-X)'^'' ,q. 

^^^'^^ {l + q)(2 + q){3 + q) (l + g)(2 + g) 1 + q ' ^^ 

In the case of the power law model, eq. |3| becomes 

x-I'dxj [l-xf-^'^X^dx, (10) 

where 

K = ^ (11) 



and L = xLm in eq. ^. Eq. |T^ has the integral solution 



N{>F)^J f{xo,8-2D)x-f'dx. (12) 



2.2. Time Dilation and Energy Shifting Factors 

In an idealized scenario of two identical bursts at different redshifts, zi and Z2, their time 
dilation and energy shifting factors, T12 and £12, are both simply equal to the ratios of their 
scale factors (neglecting the effects of energy stretching which are inherent in pulse duration 
measurements ([Fenimore &: Bloom 1995 )): 



-1 1 + ^1 ,,o\ 

n2 = ei2 = 7— — • (13) 

In practice, however, measures of the scale factor are averaged over peak flux ranges and time 
dilation and energy shifting factors are determined for pairs of these ranges. Meszaros &; Meszaros 
(1996) demonstrated that such mean values of the scale factor, averaged over a peak flux range 
Fi < F < Fu, are simple functions of the integral distribution, as modeled by eqs. ^ and ^: 

il + z)iF,F.)= ^^^^^^^^^ , (14) 

where 

N{Fi,Fu) = N{> Fu) - N{> Fi). (15) 

Consequently, time dilation and energy shifting factors between two such ranges, Fi^i < Fi < Fi^u 
and F2^i < F2 < F2_„, are given by 

_1 _ Nd+i{Fi^i,Fi^u)Nd{F2,1,F2,u) . . 

^'' '12 Nd{Fi^I,Fi,u)Nd+i{F2,i,F2,u)' ^ ' 

The effects of energy stretching are not modeled here because they are removed empirically 
from the pulse duration measurements of Norris et al. (1996a) in §3.2. The interpulse duration 
measurements of Norris et al. (1996b) and the peak energy measurements of Mallozzi et al. (1995) 
do not require such corrections. 



3. Data Analysis 
3.1. Integral Distribution 



BATSE's sensitivity becomes less than unity at peak fluxes below ~ 1 ph cm s ( [Fenimore 



et al. 1993| ). Petrosian, Lee, &: Azzam (1994) demonstrated that BATSE is additionally biased 



against short duration bursts: BATSE triggers when the mean photon count rate, defined by 

1 pt+At 

^(t) = ^/ C{t)dt (17) 



where At = 64, 256, and 1024 ms are BATSE's predefined tiniescales, exceeds the threshold count 
rate, Cum, on a particular timescale. Consequently, peak photon count rates are underestimated 
for bursts of duration T ^ At, sometimes to the point of non-detection. Peak fluxes are similarly 
underestimated. Petrosian & Lee (1996a) developed (1) a correction for BATSE's measured peak 
fluxes and (2) a non-parametric method of correcting BATSE's integral distribution. 

A burst's corrected peak flux is given by 

where F is the burst's measured peak flux and Tqq is the burst's 90% duration. Consequently, if 
Tgo ^ At, F ~ F; however if Tqq ^ At, F > F. Petrosian & Lee (1996a) demonstrated that eq. 



18| adequately corrects BATSE's measured peak fluxes (1) on the 1024 ms timescale, (2) for bursts 



of duration Tgo > 64 ms, and (3) for a variety of burst time profiles. 
BATSE's corrected integral distribution is given by 

1 (i = 1) 

I}=2(1 + ^ 



^(>^^^= nwd + i:) (^>i)' ^''^ 



where Fi > Fj+i, Fi > Fiim,i(79o), and Mi is the number of points in the associated set 
Mi = {{Fj,Fiim,jiTgo)) : Fj > Fi and ^(^^^(rgo) < F^}. The corrected threshold flux, FumiTgo), 
is the minimum value of the corrected peak flux that satisfies the trigger criterion: F > Fum, 
where 

- fF\ 

Flim — Clim -^ (20) 



and C is the measured peak photon count rate. By eq. 15, FumiTgo) is indeed a function of Tc 



tgo 



and is similarly given by 

At 

Flim{T90) = Fumil + 1^)- (21) 

-tgo 

We apply the peak flux and integral distribution corrections of Petrosian & Lee (1996a) with 
one restriction: Kouvelioutou et al. (1993), Petrosian, Lee, & Azzam (1994), and Petrosian &: 
Lee (1996a) have demonstrated that the distribution of BATSE burst durations is bimodal, with 
the division occuring at Tgo ~ 2 s. This suggests that short (Tgg < 2 s) and long (Tgg > 2 s) 
duration bursts may be drawn from separate populations. This notion is further supported by 
the tendency of short duration bursts (1) to have steeper integral distributions than long duration 



bursts ( Petrosian &: Lee 1996a ), and (2) to have lower energy shifting factors than long duration 
bursts, especially at low peak fluxes (pVIallozzi et al. 1995|) . Consequently, we exclude short 
duration bursts from our sample. 

Of the 1122 bursts in the 3B catalog, information sufficient to perform these corrections, 
subject to the above restriction, exists for 423 bursts. The corrected integral distribution is 
plotted in fig. 1. It can be seen that the corrected distribution differs significantly from the 



uncorrected distribution only at peak fluxes below F i^ 0.4 ph cm^^ s^^. For purposes of fitting, 
we truncate and normalize the integral distribution at F = 0.316 ph cni^^ s^^, which corresponds 
to a trigger efficiency of -^ g- The remaining 397 bursts are divided into eighteen bins: fifteen are 
of logarithmic length 0.1, and the brightest three are of logarithmic length 0.2. 

3.2. Time Dilation and Energy Shifting Factors 

The pulse duration time dilation factors of Norris et al. (1996a), computed using both peak 
alignment and auto-correlation statistics, are subject to energy stretching: pulse durations tend to 



be shorter at higher energies (Fenimore et al. 1995); consequently, pulse duration measurements of 



redshifted bursts are necessarily underestimated. Furthermore, Norris et al. (1996a) demonstrated 
that the unavoidable inclusion of the interpulse intervals in these analyses has a similar effect. To 
correct for these effects, Norris et al. (1996a) provided a means of calibration: they stretched and 
shifted, respectively, the time profiles and the energy spectra of the bursts of their reference bin 
by factors of 2 and 3, and from these "redshifted" bursts, they computed "observed" time dilation 
factors. For each statistic, we have fitted these "observed" time dilation factors to the "actual" 
time dilation factors of 2 and 3 with a power law which necessarily passes through the origin. 
Calibrated time dilation factors are determined from these fits and are plotted in fig. 2. 

These calibrated time dilation factors are consistent with both the interpulse duration time 
dilation factors of Norris et al. (1996b) and the energy shifting factors (long duration bursts only) 
of Mallozzi et al. (1995) (see §4), neither of which require significant energy stretching corrections. 
The interpulse duration time dilation factors were computed for various combinations of temporal 
resolutions and signal-to-noise thresholds. Norris et al. (1996b) provided error estimates for 
two such combinations, which they described as "conservative" with respect to their statistical 
significance. These time dilation factors and the energy shifting factors of Mallozzi et al. (1995) 
are additionally plotted in fig. 2. All 22 of the time dilation and energy stretching factors are fit 
to in S4. 



4. Model Fits 

Both the standard candle luminosity model and the power law luminosity distribution model 
have been x^-fitted to the corrected and binned differential distribution of fig. 1 (see §3.1) and to 
the time dilation and energy shifting factors of fig. 2 (see §3.2). Additionally, both models have 
been x^-fitted to the union of these data sets. In the case of the standard candle model, Ax^ 
confidence regions, as prescribed by Press et al. (1989), are computed on a 100^-point grid. In the 
case of the power law model, Ax^ confidence regions are computed on a 50'*-point grid and are 
projected into three two-dimensional planes. 



4.1. Standard Candle Luminosity Model 

The standard candle model consists of three parameters: h'^L^, D, and zm, where h = Hq/100. 
By eqs. ^ and ^ zm is constrained by 

->(-f(^)T-'- <"> 

where Fm = 0.201 ph cm~^ s~^ is the peak flux of BATSE's faintest burst. However, above this 
limit, Zm is independent of the data. 

The standard candle model fits both the differential distribution {xm = 18.3, v = 16) and 
the time dilation and energy shifting factors {Xm ~ 16-2, i^ = 20). The significance of the latter 
fit testifies to the consistency of the independent time dilation and energy shifting measurements. 
The Ax^ confidence regions of these fits (fig. 3), while demonstrating strong correlations between 
/i^-Lq and D, do not place bounds on either parameter. However, the latter fit places strong 
bounds on /i^Lq for reasonable values of D. 

The standard candle model additionally fits the union of these data sets {Xm = 38.2, u = 38). 
The Ax^ confidence region of this joint fit (fig. 4) places strong bounds on both h'^Lo and D: 
h^Lo = 2.3to:? X 10^^ ph s"^ and D = S.Gtoi- By eq. M, this implies that zm > Q-Ot\i, of which 



-0.7 ^ '^^ P^^ ^ '^^^"-^ ^ ~ "J-^-O.S- ^y "^H- |±f]) ''^^i^ liiipiico LiidL ^M ^ u.u_| 3, 

the implications are discussed in §5 



4.2. Po^ver La^v Luminosity Distribution Model 

The power law model consists of five parameters: h'^L, D, /?, K, and zm, where 

is the mean luminosity of the luminosity distribution, (j){L). The fifth parameter, zm, is again 
constrained by eq. ^, except with Lq -^ L^- However, unlike in the standard candle model, zm 
is not necessarily independent of the data above this limit. For purposes of fitting, we assume that 
Zm is indeed beyond what BATSE observes. The limitations of this assumption are discussed in 
§5. 

The power law model fits the differential distribution {xm = 11.2, v = 14), the time dilation 
and energy shifting factors {Xm = 13.6, u = 18), and the union of these data sets {Xm = 34.1, 
u = 36). The Ax^ confidence region of the joint fit (fig. 5) places strong bounds on D: D = 
3.7lo]5 and for h^L < 10^^ ph s~^, 3.4 ^ D ^ 3.8 to 1-a. This region is additionally divisible into 
four unique subregions (see tab. 1). Using the terminology of Hakkila et al. (1995, 1996), the 
luminosity distribution of each subregion is described as L^ dominated (independent of Lm), Lm 
dominated (independent of L^), range dominated (dependent upon both Lm and Lm), or similar 



to a standard candle {Lm ~ Lm)- For each subregion, bounds are placed on L, /?, K, and Kqq, 
where Kqq is the 90% width of the observed luminosity distribution and is given by (following the 
convention of Ulmer & Wijers (1995)) 

K,o = ^, (24) 

where Lp, the "p% luminosity" of this distribution, is defined by 

A^L<L,(> Fm) = (^) NL<L,,i> Fm). (25) 

It is important to note that others (e.g., Horack, Emslie, & Meegan (1994)) define K<jq differently: 

K - / fe^ (^"* dominated) 

^90 - I L^ (^^^ dominated) ' ^^^^ 

which results in reduced values. The former definition is applied here. 



5. Conclusions 

Assuming no evolution [D = 3), Fenimore & Bloom (1995), Nemiroff et al. (1996), and 
Horack, Mallozzi, & Koshut (1996) have demonstrated that BATSE's differential distribution is 
inconsistent with a time dilation factor of ~ 2 between the peak flux extremes of Norris et al. 
(1996a, 1996b). This has prompted suggestions that either the bursts' observed time dilation is 
largely intrinsic or that strong evolutionary effects are present in the differential distribution. 
The former explanation, however, is discredited by the degree to which the time dilation and 
energy shifting measurements are consistent. Hakkila et al. (1996), also assuming no evolution, 
have demonstrated that the differential distribution alone is incompatible with a standard candle 
luminosity. These results agree with our results for D = 2>. We additionally determine at what 
values of D that these incompatibilities disappear: D = 3.61^q 3 for the standard candle model and 
D = 3.71q'5 for the power law model. For mean luminosities h?L < 10^"^ ph s~^, evolution is even 
more tightly constrained: 3.4 ^ D <^ 3.8 (to 1-cr). 

Horack, Emshe, k Meegan (1994), Emslie & Horack (1994), Ulmer & Wijers (1995), Hakkila 
et al. (1995, 1996), and Ulmer, Wijers, & Fenimore (1995) have demonstrated that Kqq ^ 10 for 
a variety of galactic halo and cosmological models. When cosmological, these models assume no 
evolution. However, when D > 3, Kqq need not be so tightly constrained (Horack, Emslie, & 



Hartmann 1995, [Horack et al. 19961) . We find that for lO^'^ ph s"^ ^ h'^L ^ lO^'^'^ ph s-\ Kgo is 
only constrained to be less than ~ 10^ (see fig. 5). Furthermore, for h^L ^ 10^^ ph s^^, Kgo ^ 10. 
The former result is more conservative than estimates which assume no evolution. The latter is 
the result of new solutions which do not fit the data for D = 3. 



-1.5 

'1.3 

greater than that which is measured for galaxies. The power law model, under certain conditions 



In the standard candle model, the redshift of BATSE's faintest burst is 6.0_]^'3, which is much 



provides more reasonable estimates. In tab. 2, l-o" bounds are placed on the redshift of BATSE's 
faintest burst for three representative luminosities: Lio, L^q, and L90, where Lp is as defined in 
eq. |2^. (For example, L50 is the median luminosity of the observed luminosity distribution, and 
80% of the observed bursts have luminosities between Liq and -Lgo-) Defining the redshift Zp as the 
maximum redshift at which bursts of luminosity Lp can be detected, we find that 2.9 ^ 2:50 ^ 4.6 
for h^L ^ 10^^ ph s^^ and 4.2 ^ 250 ^ 9.4 otherwise. However, ziq ^ 4.2 for all mean luminosities 
and ^ 2.3 for h'^L ^ 10^^ ph s"^. If Lp ^ L90, the redshift of this burst is again quite large. 
Consequently, a mean luminosity of h^L ^ 10^^ ph s~^ coupled with a luminosity for BATSE's 
faintest burst of Lp < L50 is favored. 

In conclusion, the results presented in this paper demonstrate that when both the differential 
distribution and the time dilation and energy shifting factors are fitted to, moderate evolution is 
required if an Einstein-de Sitter cosmology, a power law spectrum of photon number index -1, 
no luminosity evolution, and in the case of the power law model, a non-observable maximum 
burst redshift are assumed. We have additionally demonstrated that under these conditions, the 
90% width of the observed luminosity distribution is not necessarily ^ 10, as appears to be the 
case if no evolution is assumed. Finally, redshift considerations indicate that if the redshifts of 
the faintest bursts are to be compatible with that which is currently known about galaxies, the 
standard candle model is unacceptable and for the power law model, a mean burst luminosity 
h^L ^ 10^^ ph cm~^ s~^ is favored. 

This work was supported in part by NASA grant NAG5-2857 and an AAS/NSF-REU grant. 
We are also grateful to E. E. Fenimore and E. D. Feigelson for useful discussions. 
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Table 1. Power Law Model Ax^ Confidence Subregions 



Subregion 



1 

2 
3 

4 



HL) 



L 



/3 



K 



Lm dominated 



K, 



90 



Lm dominated ^ Lq unbounded^ ^ 10^ ^ ]^o°-^^ 

range dominated ~ Lq ^1-5 ^ 10^ ^ 10^ 

standard candle ~ Lq unbounded ~ 1 ~ 1 

^ Lo ;^ 2.5 ^ 102-5 ^ 10 



**< 2 for cosmological values of L 
< 10^ for cosmological values of L 
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Table 2. Power Law Model Redshift of BATES's Faintest Burst'' 

Lp h^L ^ lO^'^ ph s-i h'^L ;> lO^^^ ph s^^ 

^10 1.0 ^ ziQ <; 2.3 1.2 <; zio ^ 4.2 

L50 2.9 ^ Z50 ^ 4.6 4.2 ^ Z50 ^ 9.4 

L90 5.1 ^ Z90 ^ 6.1 5.3 ^ zgo ^ 13.1 

^to l-a 
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Fig. 1. — Uncorrected and corrected integral distributions of long duration {Tqq > 2 s) BATSE 
bursts. 
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Fig. 2. — Calibrated (§3.2) pulse duration time dilation factors of Norris et al. (1996a) computed 
using peak alignment (closed circles) and auto-correlation (open circles) statistics; interpulse 
duration time dilation factors of Norris et al. (1996b) computed using temporal resolutions of 
512 ms (closed squares) and 128 ms (open squares) and signal-to-noise thresholds of 1400 counts 
s^^ (closed squares) and 2400 counts s~^ (open squares); and inverse peak energy shifting factors of 
Mallozzi et al. (1995) (closed triangles). The time dilation factors and the energy shifting factors 
have been computed using two different reference bins. 
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Fig. 3. — Ax^ confidence regions of the standard candle model fit to the differential distribution 
(A) and to the time dilation and energy shifting factors (B). Dotted lines are l-o", short dashed 
lines are 2-cr, and long dashed lines are 3-0". 
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Fig. 4. — A^^ confidence region of the standard candle model fit to the union of the differential 
distribution and the time dilation and energy shifting factors. 1-, 2-, and 3-0" are as described in 
fig. 3. 
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Fig. 5. — Projected Ax^ confidence regions of the power law model fit to the union of the differential 
distribution and the time dilation and energy shifting factors. 1-, 2-, and 3-0" are as described in 
fig. 3. Subregions 1-4 are described in tab. 1. 



